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Abstract 

The Yang-Lee zeros of the Q-state Potts model are investigated in 1, 2 and 
3 dimensions. Analytical results derived from the transfer matrix for the 
one-dimensional model reveal a systematic behavior of the locus of zeros as a 
function of Q. For 1 < Q < 2 the zeros in the complex x = exp{(3Hq) plane 
lie inside the unit circle, while for Q > 2 they lie outside the unit circle for 
finite temperature. In the special case Q = 2 the zeros lie exactly on the unit 
circle as proved by Lee and Yang. In two and three dimensions the zeros are 
calculated numerically and behave in the same way. Results are also presented 
for the critical line of the Potts model in an external field as determined from 
the zeros of the partition function in the complex temperature plane. 

I. INTRODUCTION 

The Q-state Potts model in two and three dimensions exhibits a rich variety of critical 
behavior and is very fertile ground for the analytical and numerical investigation of first- 
and second-order phase transitions. With the exception of the two-dimensional Q = 2 Potts 
(Ising) model in the absence of an external magnetic field 0, exact solutions for arbitrary Q 
are not known. However, some exact results have been established for the two-dimensional 
Q-state Potts model. For Q < 4 there is a second-order phase transition, while for Q > 4 
the transition is first order [Q. From the duality relation the critical temperature is known 
to be Tc = J/ksln^l + \/Q) [|l|]. For Q < A the critical exponents are known, while for 
Q > A the latent heat [Q, spontaneous magnetization 0, and correlation length at Tc 
are also known. 

The Q-state Potts model on a lattice G in an external magnetic field Hq is defined by 
the Hamiltonian 

nQ = JY.\l- a,)] -HqY. 5{au, q), (1) 

{ij) k 
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where J is the couphng constant, (i, j) indicates a sum over nearest-neighbor pairs, ai = 
1, Q, and g is a fixed integer between 1 and Q. The partition function of the model can 
be written as 

ZQ{x,y)= E E^q(^'^)^*V, (2) 

A/=0 E=0 

where x = e^^'', y = = e~^"^, E and M are positive integers < E < Ni, and 
< M < A^s, respectively, Nb and Ns are the number of bonds and the number of sites on 
the lattice G, and Qq{M, E) is the number of states with fixed E and fixed M. From Eq. 
(2) it is clear that ZQ{x,y) is simply a polynomial in x and y. 

By introducing the concept of the zeros of the partition function in the complex magnetic- 
field plane (Yang-Lee zeros), Yang and Lee proposed a mechanism for the occurrence of 
phase transitions in the thermodynamic limit and yielded a new insight into the unsolved 
problem of the two-dimensional Ising model in an arbitrary nonzero external magnetic field. 
It has been shown that the distribution of the zeros of a model determines its critical 
behavior. Lee and Yang ^ also formulated the celebrated circle theorem which states that 
the Yang-Lee zeros of the Ising ferromagnet lie on the unit circle. While we lack the circle 
theorem of Lee and Yang to tell us the location of the zeros for Q ^ 2, something can be 
said about their general behavior as a function of temperature. At zero temperature {y = 0) 
from Eq. (2) the partition function is 

ZQ{x,o) = J2nQ{M,o)x'' 

M 

= (Q-l)+x^^ (3) 

Therefore, the Yang-Lee zeros at T = are given by 

Xk = {Q- 1)'/^= exp[2(2A; - l)7r/Ar,], (4) 

where k = 1, ...,Ns. The zeros at T = are uniformly distributed on the circle with radius 
(Q — 1)^/^= which approaches unity in the thermodynamic limit, independent of Q. At 
infinite temperature {y = 1) Eq. (2) becomes 

Ns Nb 

Zq{x,1)= E T.^QiM,E)x^'. (5) 

M=0 E=0 

Because J2e ^q{M, E) is simply {Q — 1)^=-*^^ at T = oo, the partition function is given 
by 

ZQ(x,l) = (g-l + x)^% (6) 
and its zeros are A^^^-degenerate at x = 1 — Q, independent of lattice size. 

II. YANG-LEE ZEROS IN ONE DIMENSION 

For the one- dimensional Potts model in an external field the eigenvalues of the transfer 
matrix were found by Glumac and Uzelac 0- The two dominant eigenvalues are X± = 



2 



{A ± iB)/2a, where A ^ a{l + x) + Q - 2, B ^ -i^[a(l - x) + Q - 2]"^ + A{Q - l)x, and 
Ao = (a — is (Q — 2)-fold degenerate. 

We will assume that |A±| > Aq and verify this assumption a posteriori. The partition 
function is 



but, by the above approximation, for large N we have 

Zjv ~ A^ + A_ . 

If we define A = 2C cosip and B = 2Csmilj, where C 
X± = {C /a) exp(±i'?/'), and the partition function is 



(7) 



(8) 



Zn = 2{ — \ cos Nil) 



a — 1){Q + a — l)x^ then 



(9) 



The zeros of the partition function are then given by 

2k + I 



2N 



TT, A; = 0,1, 2,..., AT- 1. 



The location of these zeros in the complex x plane is determined by the solution of 

A = 2CcosV'fe- 
This equation is quadratic in z = ^/x with roots 
1 



a - 1){Q + a - 1) cosipk ± i\/ {a - 1){Q + a - 1) sin^ Vfe + <5 - 1 



It is easily verified that 



Q + a-2 



(10) 



(11) 



(12) 



(13) 



and we see that all the zeros lie on a circle in the complex x plane. The argument of Xk is 
given by 



cos — = 



(a-l)(g + a-l) 



2 \ a(Q + a-2) 



cos V'fe. 



(14) 



Before we analyze these results we must verify that in fact |A±| > Aq in the region of 
these zeros. We find 



I A 



± 



Ao \ 



(g + g- l)(g + a-2) 
a{a — 1) 



(15) 



which is indeed greater than unity for g > 1. 
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If we return to Eq. (13) we can discern a remarkably systematic behavior of the Yang-Lee 
zeros with Q. For 1<Q<2, |xfc|<l and the zeros he inside the unit circle. As a — cxd, 
the zeros approach the unit circle from within, as we observed in Eq. (4), and for a = 1, 
cos(6'fc/2) = 0, so 6*^ = vr and all the zeros lie at 1 — Q, as predicted in Eq. (6). For Q > 2, 
\xk\ > 1 and the zeros lie outside the unit circle and approach the unit circle as a ^ oo from 

outside. In the special case Q = 2 we of course find that | X/j;| 1, clS proved by Lee and 

Yang 0. The edge singularity in the thermodynamic limit is given by 



Oq = 2 cos" 



(a-m + a-l)^^^ (16) 



\ a{Q + a-2) 



from which we conclude that no transition occurs for any T > 0. 

Finally, we can use these exact results to study finite size effects on the distribution of 
zeros. For finite N we must find the zeros of the full partition function as given in Eq. (7). 
We find that finite size effects are very small in one dimension, but the largest deviation 
from the limiting locus occurs for arg(x) = vr. 



III. YANG-LEE ZEROS IN TWO AND THREE DIMENSIONS 

We have calculated exact integer values for Qq{M,E) on L x L square lattice for 3 < 
L < 8 and 3'^ simple-cubic lattice {Q = 3) using the restricted microcanonical transfer 



matrix For square lattices up to L = 12 {Q = 3) the zeros are calculated using a 

semi- canonical variation of the transfer matrix. Fig. la shows the Yang-Lee zeros for the 
two-dimensional three-state Potts model in the complex x plane at the critical temperature 
He = 1/(1 + v^) = 0.366... for L = 4 and L = 10 with cylindrical boundary conditions. 
Note that just as in the one-dimensional model, the zeros of the three-state Potts model 
lie close to, but not on, the unit circle. The zero farthest from the unit circle is in the 
neighborhood of arg(x) = vr, while the zero closest to the positive real axis lies closest to 
the unit circle. The behavior of the zeros with the size of the lattice also follows that of 
the one- dimensional model; the zeros for L = 10 lie on a locus interior to that for L = 4, 
and we find similar behavior for larger values of Q. In the thermodynamic limit the locus of 
zeros cuts the real x axis at the point x = 1 |[ll[] corresponding to Hg = 0, as described by 



Yang and Lee 0. Fig. lb shows the zeros for the two-dimensional three-state Potts model 
at several temperatures with cylindrical boundary conditions. At y = 0.5yc the zeros are 
nearly uniformly distributed in argument and close to the unit circle. As the temperature 
is increased the edge singularity moves away from the real axis and the zeros detach from 
the unit circle. Finally, as y approaches unity, the zeros converge on the point x = —2. As 
predicted by Eqs. (4) and (6), for periodic boundary conditions |T2| and self-dual boundary 



conditions [|TT| we observe the same behaviors as those in Fig. 1 for cylindrical boundary 
conditions. 

The exact nature of the locus of zeros for the Q-state Potts model in two dimensions 
is unknown, with the exception of Q = 2. It is clear that the Yang-Lee zeros of the two- 
dimensional Q-state Potts model do not lie on the unit circle for Q > 2 for any value of 
y and any finite value of L. Since the zero in the neighborhood of arg(a;) = vr is always 
the farthest from the unit circle, if this zero can be shown to approach |a;(7r)| = 1 in the 
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limit L — > oo, all the zeros should lie on the unit circle in this limit. In Fig. 2 we show 



values for |x(7r)| extrapolated to infinite size using the Bulirsch-Stoer (BST) algorithm |13 
for 3<Q<8a.ty = 0.5?/c and y = yc- The error bars are twice the difference between the 
(n — 1, 1) and (n — 1,2) approximants. From these results it is clear that while the locus 
of zeros lies close to the unit circle at y = yc, it does not coincide with it except at the 
critical point a; = 1. While the numerical evidence presented here suggests that the locus of 
zeros for the three-state model is not the unit circle, and the exact result in one-dimension 
provides an example where the locus varies continuously with Q, the nature of the locus 
remains an open and fascinating question. 

Fig. 3 shows the BST estimate of the modulus of the locus of zeros as a function of angle 
for the two-dimensional three-state Potts model at y = 0.5?/c, y = Vc and y = 1.2?/c- To 
calculate the extrapolated values for each angle, 9, we selected the zero whose arguments 
were closest to 6 for lattices of size 3 < L < 12 for 6 = 0.0, 0.5, 2.5 and vr. The BST 
algorithm was then used to extrapolate these values for finite lattices to infinite size. The 
large variation in the size of the error bars is due to the fact that for a given 6 there may 
be no zero close to 9 for the smaller lattices. In Fig. ?> at y = y^ the first four angles are 
shifted slightly from the original values {9 = 0.0, 0.5, 1.0 and 1.5) to be distinguished from 
the results at y = 0.5j/c. For y = 0.5?/c and y = yc the first zeros definitely lie on the point 
r{9 = 0) = 1 in the thermodynamic limit. However, for y = 1.2yc the BST estimates of 
the modulus and angle of the first zero are 1.054(2) and 0.09(6). Therefore, at y = 1.2?/^ 
the locus of zeros does not cut the positive real axis in the thermodynamic limit, consistent 
with the absence of a physical singularity for y > yc- 

Fig. 4 shows the Yang-Lee zeros for the three-dimensional three-state Potts model at 
several temperatures on the 3^ simple-cubic lattices with periodic boundary conditions in x 
and y directions and free boundary conditions in z direction. The behaviors in Fig. 4 for the 
three-dimensional model are the same as those in Fig. lb for the two-dimensional model. 



IV. FISHER ZEROS IN AN EXTERNAL FIELD 

Fisher emphasized that the partition function zeros in the complex temperature 
plane (Fisher zeros) are also very useful in understanding phase transitions. In particular, 
using the Fisher zeros both the ferromagnetic phase and the antiferromagnetic phase can be 
considered at the same time. From the exact solutions of the square lattice Ising model it 
has been shown |T^ that in the absence of an external magnetic field the Fisher zeros lie on 



two circles in the thermodynamic limit. Recently the locus of the Fisher zeros of the Q-state 
Potts model in the absence of an external magnetic field has been studied extensively |]T3|,|TB[. 
It has been shown [|16[ that for self-dual boundary conditions near the ferromagnetic critical 
point yc = 1/(1 + VQ) the Fisher zeros of the Potts model on a finite square lattice lie on 
the circle with center —1/{Q — 1) and radius \fQI[Q — 1) in the complex y plane, while the 
antiferromagnetic circle of the Ising model completely disappears for Q > 2. However, the 
properties of the Fisher zeros for Q > 2 in an external field are not known. 

In the limit Hq — oo {x 0) the partition function of the Q-state Potts model 
becomes 
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ZQ=J2QQ{M = 0,E)y^, (17) 

E=0 

where Qq{M = 0,E) is the same as the number of states Qq^i{E) of the {Q — l)-state 
Potts model in the absence of an external magnetic field. As x decreases from 1 to 0, the 
Q-state Potts model is transformed into the {Q — l)-state Potts model in zero external field 
||T7| . For an external field Hg < 0, one of the Potts states is supressed relative to the others. 
The symmetry of the Hamiltonian is that of the {Q — l)-state Potts model in zero external 
field, so that we expect to see cross-over from the Q-state critical point to the {Q — l)-state 
critical point as —Hg is increased. 

We have studied the field dependence of the critical point for < s < 1 through the 
Fisher zero closest to the real axis, yi{x, L), for the two-dimensional three-state Potts model. 
For a given applied field yi approaches the critical point yi{L) ydx) in the limit L — > cxo, 
and the thermal exponent yt{L) defined as ||TO| , p^ 



.n- Mlmbi(L + l)]/Im[yi(L)]} 

- ln[(L + 1)/L] ^^^^ 

will approach the critical exponent yt{x). Table I shows values for y^x) extrapolated from 
calculations of yi{x,L) on L x L lattices for 3 < L < 8 using the BST algorithm. The critical 
points for x = 1 (three-state) and a; = (two-state) Potts models are known exactly and 
are included in Table I for comparison. Note that the imaginary parts of ?/c(BST) are all 
consistent with zero. We have also calculated the thermal exponent, yt, applying the BST 
algorithm to the values given by Eq. (18), and these results are also presented in Table I. 
For X = 1 we find yt very close to the known value yt = 6/5 for the three-state model, but 
for X as large as 0.5 we obtain yt = 1, the value of the thermal exponent for the two-state 
(Ising) model. 

Fig. 5 shows the critical line of the two-dimensional three-state Potts ferromagnet for 
Hq < 0. In Fig. 5 the upper line is the critical temperature of the two-state model, 
Tc{Q = 2) = l/ln(l -|- V2), and the lower line is the critical temperature for the three- 
state model, Tc{Q = 3) = l/ln(l -|- v^). The critical line for small —Hg is given by 
T - T^{Q = 3) ~ {-Hq)y'/y\ where = 6/5 and yn = 28/15 for the three-state Potts 
model. 



V. CONCLUSION 

Following the work of Glumac and Uzelac ||^ we have studied the locus of Yang-Lee 
zeros in the complex x plane for the Q-state Potts model in one-dimension and find that 
for Q > 1 the zeros lie on a circle whose radius varies continuously with temperature and 
Q. For 1 < Q < 2 and any finite temperature, the zeros lie inside the unit circle, while for 
Q > 2 they lie outside the unit circle. In the special case Q = 2 (Ising model) the zeros lie 
exactly on the unit circle, as first proved by Lee and Yang [0. 



In two and three dimensions we have used the microcanonical transfer matrix [|T0| to 
find the Yang-Lee zeros for finite lattices. The general trends observed in one dimension 
are repeated in two and three dimensions. In two dimensions with 3 < Q < 8 the zeros 
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lie outside the unit circle, and finite-size scaling suggests that the locus of zeros in the 
thermodynamic limit touches the unit circle only at the critical point, x = 1. In three 
dimensions with Q = 3 we have calculated the zeros on a single 3x3x3 lattice, and there 
again the zeros lie outside the unit circle. While the exact form of the locus of zeros remains 
an open and important question, our results suggest a universal behavior which incorporates 
the Lee- Yang circle theorem. For 1 < Q < 2 the locus of zeros lies within the unit circle 
while for Q > 2 it lies outside the unit circle and in each case touches the unit circle at the 
critical point x = 1. 

We have also studied the Fisher (complex temperature) zeros of the two-dimensional 
three-state Potts model for lattices of size 3 < L < 8 in the presence of an external field. 

The field reduces the symmetry of the Hamiltonian to that of the {Q — l)-statc model, and 
by studying the edge singularity we are able to determine the critical point and temperature 
exponent as a function of the field. Cross-over from the 3-state to the 2-state universality 
class is apparent. A more detailed analysis using larger lattices and allowing for confluent 
singularities is in progress. 
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TABLES 

TABLE I. The critical temperature yc and the critical exponent yt of the two-dimensional 
three-state Potts model for < x < 1. 



X Vc (BST) yc (exact) yt (EST) yt (exact) 

0.414(3) +0.0002(4)i 0.414213... 1.001(2) 1 
0.001 0.414(3) + 0.0002(4)i 1.001(2) 

0.05 0.413(5) + 0.0002(3)z 1.0009(6) 

0.5 0.400(2) + 0.000(2)z 0.982(21) 

1 0.366(2) + 0.0002(5)z 0.366025... 1.195(3) | 
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FIG. 1. Zeros of the two-dimensional three-state Potts model in the complex x plane with 
cylindrical boundary conditions (a) at y = yc ioi L = 4 and L = 10 and (b) for several values of y 
(L = 6) 
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FIG. 2. Modulus of the zero at 9 = tt extrapolated to infinite size for 3 < Q < 8 at y = 0.5?/c 
and y = Vc with cylindrical boundary conditions. 
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CD 




FIG. 3. Modulus of the locus of zeros as a function of angle for the two-dimensional three-state 
Potts model at y = 0-5yc, Vc, and 1.2yc with cylindrical boundary conditions. The slight horizontal 
off-set for data for y = yc'is for clarity only. However, the off-set of the edge singularity for y = 1.2j/c 
from ^ = is real. 
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FIG. 5. Critical temperatures of the two-dimensional three-state Potts ferromagnet as a func- 
tion of the magnetic field. Hq is in unit of J and T is in unit of J/ks- The upper dotted line is the 
Ising transition temperature in the limit Hq ^ — oo, while the lower dotted line shows the critical 
temperature of the three-state Potts model for Hq = 0. 
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